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Abstract 

We present models of flavorful supersymmetry in higher dimensions. The Higgs fields 
and the supersymmetry breaking field are localized in the same place in the extra di- 
mension^). The Yukawa couplings and operators generating the supersymmetry breaking 
parameters then receive the same suppression factors from the wavefunction profiles of the 
matter fields, leading to a specific correlation between these two classes of interactions. 
The resulting phenomenology is very rich, while stringent experimental constraints from 
the low-energy flavor and CP violating processes can all be satisfied. We construct both 
unified and non-unified models in this framework, which can be either strongly or weakly 
coupled at the cutoff scale. We analyze one version in detail, a strongly coupled unified 
model, which addresses various issues of supersymmetric grand unification. The models 
presented here provide an explicit example in which the supersymmetry breaking spectrum 
can be a direct window into the physics of flavor at a very high energy scale. 



1 Introduction 



One of the longstanding puzzles of the standard model is the distinct pattern of masses and 
mixings of the quarks and leptons. While supersymmetry addresses many of the other mysteries 
of the standard model, including the instability of the electroweak scale and the lack of a dark 
matter candidate, it is not clear if and how supersymmetry helps us understand the flavor puzzle 
of the standard model at a deeper level. Recently, it has been pointed out that the supersymme- 
try breaking parameters can exhibit nontrivial flavor structure, and that measurement of these 
parameters at the LHC can give insight into the flavor sector of the standard model [1, 2]. In 
particular, it has been shown in Ref. [2] that the class of models called flavorful supersymme- 
try, in which the supersymmetry breaking parameters receive similar suppressions to those of 
the Yukawa couplings, can evade all the current experimental bounds and have very distinct 
signatures at the LHC. In this paper we present explicit models of flavorful supersymmetry. 

In this paper we construct models in higher dimensional spacetime where supersymmetry 
breaking and the Higgs fields reside in the same location in the extra dimension(s). This provides 
a simple way to realize the necessary correlation between the structures of the supersymmetry 
breaking parameters and the Yukawa couplings [3, 2]. To preserve the successful prediction for 
super symmetric gauge coupling unification, we take the size of the extra dimension(s) to be of 
order the unification scale. The hierarchical structure for the Yukawa couplings is generated by 
wavefunction overlaps of the matter and Higgs fields [4], and the correlation between flavor and 
supersymmetry breaking is obtained by relating the location of the Higgs and supersymmetry 
breaking fields in the extra dimension(s). Models along similar lines were considered previously 
in Ref. [5], where flavor violation in the supersymmetry breaking masses is induced by finite 
gauge loop corrections across the bulk. Here we consider models in which matter fields interact 
directly with the supersymmetry breaking field, giving the simplest scaling for flavorful effects 
in the supersymmetry breaking parameters. 1 

While not necessary, the extra dimension(s) with size of order the unification scale can also 
be used to address various issues of supersymmetric grand unified theories. Grand unification in 
higher dimensions provides an elegant framework for constructing a simple and realistic model 
of unification [8, 9]. It naturally achieves doublet-triplet splitting in the Higgs sector and sup- 
presses dangerous proton decay operators, while preserving successful gauge coupling unification. 
Realistic quark and lepton masses and mixings are also accommodated by placing matter fields 
in the bulk of higher dimensional spacetime [9, 10, 11]. We thus first construct a grand unified 
model of flavorful supersymmetry which can successfully address these issues. In this model 

1 Flavor violation in higher dimensional supersymmetric models was also discussed in different contexts, see 
[6, 7]. 
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we also adopt the assumption of strong coupling at the cutoff scale motivated by the simplest 
understanding of gauge coupling unification in higher dimensions [12, 13], although this is not a 
necessity to realize flavorful supersymmetry. 

There are a variety of ways to incorporate supersymmetry breaking in the present setup. 
An important constraint on the flavorful supersymmetry framework is that superpotential op- 
erators leading to the supersymmetry breaking scalar trilinear interactions must be somewhat 
suppressed, unless the superparticles are relatively heavy. While it is possible that this suppres- 
sion arises accidentally or from physics above the cutoff scale, we mainly consider the case where 
the suppression is due to a symmetry under which the supersymmetry breaking field is charged. 
This symmetry can also be responsible for a complete solution to the // problem, the problem of 
the supersymmetric Higgs mass term (the \x term) being of order the weak scale and not some 
large mass scale. This leads to a scenario similar to the one discussed in Refs. [14, 15], in which 
the \i term arises from a cutoff suppressed operator [16] while the gaugino and sfermion masses 
are generated by gauge mediation [17, 18]. The present setup, however, also leads to flavor vio- 
lating squark and slepton masses that are correlated with the Yukawa couplings, characterizing 
flavorful supersymmetry. 

We stress that only the extra dimension(s) and the field configuration therein are essential 
for a realization of flavorful supersymmetry. All the other ingredients, including grand unifi- 
cation, strong coupling, and the particular way of mediating supersymmetry breaking, are not 
important. While the model described above provides an explicit example of flavorful supersym- 
metry in which many of the issues of supersymmetric unification are addressed in a relatively 
simple setup, it is straightforward to eliminate some of the ingredients or to extend the model 
to accommodate more elaborate structures. In particular, we explicitly discuss a construction in 
which the theory is weakly coupled at the cutoff scale, which can be straightforwardly applied 
to models with various spacetime dimensions or gauge groups. 

The organization of the paper is as follows. In the next section we present a unified model of 
flavorful supersymmetry with the assumption that the theory is strongly coupled at the cutoff 
scale. We explain how the relevant correlation between the Yukawa couplings and supersymmetry 
breaking parameters is obtained. Phenomenology of the model is studied in section 3, including 
constraints from low-energy processes, the superparticle spectrum, and experimental signatures. 
In section 4 we construct a model in warped space, which allows us to obtain a picture of 
realizing flavorful supersymmetry in a 4D setup, through the AdS/CFT correspondence. In 
section 5 we present a weakly coupled, non-unified model, which does not possess a symmetry 
under which the supersymmetry breaking field is charged. Extensions to larger gauge groups or 
higher dimensions are also discussed. Finally, conclusions are given in section 6. 
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2 Model 

In this section we present a unified, strongly coupled model. We adopt the simplest setup, SU (5) 
in 5D, to illustrate the basic idea. Extensions to other cases such as larger gauge groups and/or 
higher dimensions are straightforward. It is also easy to reduce the model to a non-unified model 
in which the gauge group in 5D is the standard model SU(3)c x SU{2)l x U(1)y- 

2.1 SU(5) grand unification in 5D 

We consider a supersymmetric SU (5) gauge theory in 5D flat spacetime with the extra dimension 
compactified on an S 1 /Z 2 orbifold: < y < nR, where y represents the coordinate of the extra 
dimension [8, 9]. Under 4D N = 1 supersymmetry, the 5D gauge supermultiplet is decomposed 
into a vector superfield V(A jX , A) and a chiral superfield E(<7 + iA$, A'), where both V and E are 
in the adjoint representation of SU(5). We impose the following boundary conditions on these 
fields: 
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where + and — represent Neumann and Dirichlet boundary conditions, and the first and second 
signs in parentheses represent boundary conditions at y = and y = nR, respectively. This 
reduces the gauge symmetry at y — nR to £77(3) x £77(2) x U(l), which we identify with the 
standard model gauge group SU(3)c x SU(2)l x U(1)y (321). The zero-mode sector contains 
only the 321 component of V, V 321 , which is identified with the gauge multiplet of the minimal 
supersymmetric standard model (MSSM). 

The Higgs fields are introduced in the bulk as two hypermultiplets transforming as the 
fundamental representation of ££7(5). Using notation where a hypermultiplet is represented by 
two 4D N = 1 chiral superfields $(0,^) and $ C (0 C , ip c ) with opposite gauge transformation 
properties, our two Higgs hypermultiplets can be written as {H, H c } and {H, H c }, where H and 
H c transform as 5 and H and H c transform as 5* under SU(5). The boundary conditions are 
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given by 

if (5) = H T (3,l)% ] ®H D (l,2) ( j 2 + \ (3) 

ff c (5*) = ff^(3M)S7 3 +) ©ff| ) (i,2) ( : 1 7 2 ) , (4) 

for {H,H C }, and similarly for {H,H C }. Here, the right- hand- side shows the decomposition of 
H and H° into representations of 321 (with U(1) Y normalized conventionally), together with 
the boundary conditions imposed on each component. The zero modes consist of the SU(2) L - 
doublet components of H and H , H D and H D , which are identified with the two Higgs doublets 
of the MSSM, H u and H d . 

Matter fields are also introduced in the bulk. To have a complete generation, we introduce 
three hypermultiplets transforming as 10, {T, T c }, {T',T' C } and {T",T" C }, two transforming 
as 5*, {F, F c } and {F',F' C }, and one transforming as 1, {0,O c }, for each generation. The 



boundary conditions are given by 

T(10) = T Q (3,2)S; 6 +) ©T l/ (3M) (+ 2/ - ) ©T s (l,l)S + '" ) , (5) 

T'(10) = T4(3,2)S;- ) ©T/ / (3%l)L + 2/ + 3 ) ©^(l,l)S + '- ) , (6) 

T"(10) = T^(3,2)S; 6 - ) ©T/;(3*,l)L + 2 73 ) ©^(l ) l)S + ' +) , (7) 

F(5*) = F D (3\1)[% +) ®F L (1,2)%\ (8) 

F'(5*) = FUS\1)^ 3 ~ } ®F^1,2)%\ (9) 

0(1) = O^l,!)^. (10) 



The boundary conditions for the conjugated fields are given by + <-> — , as in Eqs. (3, 4). With 
these boundary conditions, the zero modes arise only from Tq, T{j, T£, Fd, F' l and On, which 
we identify with a single generation of quark and lepton superfields of the MSSM (together with 
the right-handed neutrino), Q, U, E, D, L and N 2 

There are two important scales in the theory: the cutoff scale M* and the compactification 
scale 1/R. We take the ratio of these scales to be irRM* lQir 2 /g 2 C xs O(10 - 100), where g is 
the 4D gauge coupling at the unification scale, g = 0(1), and C ~ 5 is the group theoretical factor 
for SU(5). This makes the theory strongly coupled at M*, suppressing incalculable threshold 
corrections to gauge coupling unification [12, 13]. 3 Motivated by successful gauge coupling 

2 It is possible to extract both U and E from a single hypermultiplet {T',T' C } by adopting the boundary 
conditions T'(10) = T^{3,2)[+f ] T^(3*, 1)^/3 ® t e0-> in which case we do not introduce the 

hypermultiplet {T",T" C }. In fact, this is what we obtain if we naively apply the orbifolding procedure to the 
matter hypermultiplets. The model also works in this case, with the extra constraint of M\j i = (see 
section 2.2) and qq = ql (see section 2.3). 

3 Our estimate on the strong coupling scale is conservative. It is possible that M*i? can be larger by a factor 
of w 7r, but it docs not affect our results. 
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unification at about 10 16 GeV in supersymmetric models, we take the cutoff scale and the scale 
of the extra dimension to be 



M* « 10 17 GeV, 



1/nR « 10 15 GeV. 



(11) 



More detailed discussions on gauge coupling unification will be given in section 3.4. 
2.2 Quark and lepton masses and mixings 

With the boundary conditions given in the previous subsection, the matter content of the theory 
below 1/R reduces to that of the MSSM and right-handed neutrinos: V 321 , H u , Hd, Qi, Ui, Di, 
Li, Ei and N iy where i = 1,2,3 is the generation index. The Yukawa couplings for the quarks 
and leptons are introduced on the y — and y = tcR branes. The sizes of the 4D Yukawa 
couplings are then determined by the wavefunction values of the matter and Higgs fields on 
these branes. This can be used to generate the observed hierarchy of quark and lepton masses 
and mixings [4, 5, 11]. Here we consider particular configurations of these fields, relevant to our 
framework. 

A nontrivial wavefunction profile for a zero mode can be generated by a bulk mass term. A 
bulk hypermultiplet {$, $ c } can generally have a mass term in the bulk, which is written as 



in the basis where the kinetic term is given by S kin = J d A x J dy [J d 4 8 ($ t <l>+$ c $ ct )+{ / d 2 6 § c d y &+ 
h.c.}] [19]. The wavefunction of a zero mode arising from $ is proportional to e~ Mi ' y , so that it 
is localized to the y = (y = nR) brane for M$ > (< 0), and flat for M$ = 0. (The $ c case is 
the same with M$ — > — M$.) In the present model, we have a bulk mass for each of the Higgs 
and matter hypermultiplets. For clarity of notation, we specify these masses by the subscript 
representing the corresponding zero mode: M Hu , M Hi) Mq v M Uv M d . : M l . : M e . and M N .. 

We mainly consider the case that the two Higgs doublets H u and H d are strongly localized 
to the y = ttR brane: 




(12) 



M Hu , M Hd < --. 
The relevant Yukawa couplings are then those on the y = ttR brane 



(13) 



S = 




(14) 
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where the sizes of the couplings are naturally given by (\ u )ij, (^d)ij, (^e)ij, ~ 47r/M*^ 2 

using naive dimensional analysis [20, 12]. This leads to the low-energy 4D Yukawa couplings 

W = (yu)ijQiUjH u + (yd)ijQiDjH d + (y e ) ij L i E j H d + (y v )ijLiNjH v , (15) 

with 

(y u )ij « 47re Qi e t 7 j , « 47re Qi e Dj , (y e )y « 47re Li e Bj , (y^ « 4ne L .e N ., (16) 

where the factors e$ ($ = Qi,Ui, D iy Li, E iy iVj) are given by 



2M 4 



2M* e -xflM $ for £ ! 



e$ = J ^rr- e^ RM * ~ <| for (vr^M*! < 1 . (17) 

v 1 ; [ Vt for * - 1 

Realistic Yukawa couplings are obtained by localizing lighter generations more towards the 
y = brane so that their wavefunction overlaps with the Higgs fields are more suppressed. For 
example, we can take 

eQ 2 ~y~^e, £u 2 ~y~^, e D2 ^y-h, e L2 & y~h, e E . 2 ~y~^e, (18) 
e Q 3 ~^"i £u 3 ~y~K e Ds ^y-h, e L ,^y-h, e Es ~ 

tan/3 = W ^e-i, (19) 

where e ~ O(0. 1) and y ~ 47r ps 1/e, to reproduce the gross structure of the observed quark 
and lepton masses and mixings. The suppression factors of Eq. (18) are obtained by taking bulk 
masses 

1 5-1 15 

M Q 3 ,u 3 ,e 3 ~ M Q2:U2tDi:LhE2 « — — , M Ql:UuEl w — . (20) 

Small neutrino masses are obtained through the seesaw mechanism by introducing Majorana 
masses for the right-handed neutrinos on the y = nR brane 



and 



S = Jd 4 x J^dy 5(y - ttR) J d 2 9 O m O Nj + h.c. (21) 

The values of are then not relevant to the low-energy masses and mixings (unless iVj's are 
localized to the y = brane extremely strongly), since they cancel out in the expression for the 
light neutrino masses. 

The localization of various fields in the extra dimension with the bulk masses of Eqs. (13, 20) 
is depicted schematically in Fig. 1. The quark and lepton masses and mixings are given by 
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Figure 1: A schematic depiction of the localization for various fields. Here, X represents the 
super symmetry breaking field (see section 2.3). 
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(22) 



(23) 



where 0(1) factors are omitted from each element, and Vckm and Vmns are the quark and lepton 
mixing matrices, respectively. This reproduces the gross structure of the observed quark and 
lepton masses and mixings [21]. 

The matter configuration considered here can be extended easily to account for the more 
detailed pattern of the observed masses and mixings. For example, we can localize L x slightly 
more towards the y = brane to explain the smallness of the e3 element of Vmns, which is 
experimentally smaller than about 0.2. The other elements of Vckm and Vmns, as well as the 
mass eigenvalues, can also be better fitted by choosing the bulk masses more carefully. Here we 
simply adopt Eq. (20) (and its variations, discussed in section 3.1) for the purpose of illustrating 
the general idea. 

There are also variations on the location of the Higgs fields. For example, we can localize 
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the two Higgs doublets on the y — brane, instead of the y = irR brane: M Hu ,M Hd ^ 1/nR. 
In this case, the localization should not be very strong so that their colored-triplet partners, 



location of the matter fields can simply be flipped with respect to y = ttR/2: M$ — > — M$ for 
$ = Qi,Ui, Di, Li, Ei, JVj. Another possibility is to (slightly) delocalize H u and/or H d from the 
brane. In this paper, we focus on the case of Eq. (13), where H u and H d are strongly localized 
to the y = tiR brane. 

2.3 ii term, U(1)h, and flavorful super symmetry 

In order to have a complete solution to the doublet-triplet splitting problem, a possible large 
mass term for the Higgs doublets on the y = tcR brane, 5(y — irR) J d 2 9 H D H D) must be forbidden 
by some symmetry. Moreover, to understand the weak scale size of the mass term (/z term) for 
the Higgs doublets, the breaking of this symmetry must be associated with supersymmetry 
breaking. One possibility to implement this idea is to consider a U(1)r symmetry under which 
the two Higgs doublets are neutral [9] . Here we consider the case that the symmetry is a non-i? 
symmetry. 

We consider that the bare [i term, / d 2 9 H u H d , is forbidden, but the effective \x term is 
generated by the operator Jd A 6 X* H u H d through supersymmetry breaking, where X is a super- 
symmetry breaking field [16]. We then find that the relevant symmetry is U(l) (a Peccei-Quinn 
symmetry) whose charge assignment can be taken, without loss of generality, as 



QMq), Ui(-l-q Q ), A(-l-5o), HQl), Ei(-l-q L ), Ni(-l-q L ), (24) 



where qQ and qi are real numbers, and we have assumed that the Yukawa couplings are invariant 
under the symmetry. In the context of the 5D theory, this assignment can be implemented by 
considering U (1) charges for a hypermultiplet {$, $ c } ($ = H, H, T iy T/, T", F iy Oi) such that 
the charge of $ follows that of the corresponding zero mode, while the charge of $ c is the 
opposite to that of This U(l) symmetry commutes with 5D supersymmetry. The X field 
is introduced on the y = irR brane, either as a brane field or a bulk field whose zero mode is 
strongly localized to the y = irR brane by a bulk mass term Mx «C — 1/R (see Fig. 1). 

The U(l) symmetry of Eqs. (24, 25), which we call U(1)h, has several immediate virtues. 
First of all, the most general interactions between the Higgs and X fields, located on the y = ttR 
brane, leads (up to the quadratic order in X) to the following interactions in 4D: 



whose masses are given by w 2M Hu e 



and 2M Hd e- nRMn d, do not become too light. The 



H u (l), H d (l), X(2), 



(25) 




X^H u H d + h.c.) + ^XH\H U + ^XH\H d , 



(26) 
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where we have used naive dimensional analysis to estimate the sizes of various coefficients, and 
omitted an 0(1) factor in each term. The mass scale A is defined by 

A = ^ » 10 16 GeV, (27) 

47T 

where we have used Eq. (11). After supersymmetry is broken by the F-term vacuum expectation 
value (VEV), F X) of the X field (see the next subsection), these interactions lead to the fx term 
and soft supersymmetry breaking masses for the Higgs fields of order F x /A at the scale M*: 

Fx ( Fx \ 2 

(Note that 0(1) coefficients are omitted in these equations, so that the ratio of m 2 ^ to fn 2 Hd , 
for example, can be an arbitrary 0(1) number.) An important point here is that the operator 
C Jd 4 6 (X*i X H u H d / A 2 + h.c.) is prohibited by U(1) H , so that the holomorphic supersymmetry 
breaking mass-squared for the Higgs doublets (Bfx term) is not generated at order (Fx/A) 2 at tree 
level. 4 The low-energy value of the B\i term is then generated by contributions from the gaugino 
masses through renormalization group evolution. This is crucial to avoid the supersymmetric 
CP problem, since for weak scale superparticle masses an arbitrary relative phase between the 
fi and Bjj, terms leads to an unacceptably large electric dipole moment for the electron. 

Another important implication of U{1)h is that possible y = nR brane operators 5(y — 
nR) Jd 2 9 (XT Q T^H D +XT Qi F Dj H D +XF' Li T'^H D +XF' Li O Nj H D )+\i. C ., which reduce in 4D to 
/ d 2 9 (XQiUjH u +XQiDjH d +XLiEjH d +XLiNjH u )+h.c., are forbidden. If these operators were 
present, they would lead to supersymmetry breaking scalar trilinear interactions (A terms) of 
order (a u )ij ^ 47reQ i e [ /.(Fx/A), (a d )ij ^ A-Ke Qi e Dj (F x / 'A), (a e )^ ^ A-Kt Lt t Ej {F x / A) and (a v )ij ^ 
AneL^NjiFx/A), which are not necessarily proportional to the corresponding Yukawa matrices 
in flavor space. Here, (a/)ij (/ = u, d, e, v) are defined by £ so ft = —(a u )ijqiUjh u — (ad)ijqidjh d — 
{a^ijliCjhd — (a^ijlihjhu + h.c. While these terms are suppressed by e factors, they still provide 
sizable contributions to low-energy flavor violating processes, because an A-term insertion flips 
the chirality of the sfermion and thus eliminates one factor of the Yukawa coupling from an 
amplitude. We then find that with 0(1) coefficients, the rate for /i — > is expected to be 
larger than the experimental upper bound by a couple of orders of magnitude for weak scale 
superparticle masses [2, 3]. This problem does not arise in the present model. 

The interactions between the matter and X fields relevant to soft supersymmetry breaking pa- 
rameters take the form 5(y-nR) Jd 4 9 {X^ XTq\Tq .+X^ XT^T^ - +X^ X F D \F D j+X^ X F' I } i F' L - + 

4 There are contributions to the B/j, term of order F^-(X) / 'A 3 and Fxm 3 / 2 /A 7 where m 3 / 2 is the gravitino 
mass, coming from operators C w Jd 4 9 (X^ 2 XH u Hd/A 3 + h.c.) and the supergravity effects of the first term 
of Eq. (26), respectively. These contributions arc, however, negligibly small, since (X)/A w A/Mpi <C 1 and 
m 3/2 w Fx/Mpi < F x /A, where M P i ~ 2 x 10 18 GeV is the reduced Planck scale (see section 2.4). 
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X^XT^T'^ + X*XO N \0 Nj ), which reduce in 4D to 

/^EE^*!*;. ( 29 ) 

where $ = Q, U, D, L, E, N. This leads to the following supersymmetry breaking squared masses 
for the squarks and sleptons at the scale M*: 

(F \^ (F \^ /.f 1 \^ 

-^■J > ( m l)ij ~ (xj ' ^ m ^' ~ €Di€Dj vxj ' ^ 30 ^ 

(F \^ f F \ 2 

Xj , ~ e^e^ (^J , (31) 

where we have omitted supersymmetry breaking masses for the right-handed sneutrinos, which 
are not relevant for low-energy phenomenology. Through Eq. (16), these masses are related 
to the Yukawa couplings — lighter generation scalars receive only small contributions, while 
heavier generation scalars can receive sizable ones. This is exactly the pattern needed to re- 
alize the flavorful supersymmetry scenario, which arises here from the fact that the Higgs and 
supersymmetry breaking fields reside in the same location in the extra dimension. 

As shown in Ref. [2], the existence of flavor non-universal contributions of Eqs. (30, 31) 
does not contradict the low-energy data on flavor or CP violating processes for wide parameter 
regions. Since the masses of Eqs. (30, 31) are highly flavor non-universal, they cannot be the 
dominant contribution to the soft masses (except possibly for some of the third generation 
sfermions), and we need an extra flavor universal contribution as well as the gaugino masses. 
These are generated in the present model by gauge mediation, as discussed in the next subsection. 

Finally, the U(1)h symmetry forbids any superpotential term involving only the X field. 
Since breaking supersymmetry requires a linear X term in the superpotential, this implies that 
supersymmetry is not broken unless U(1)h is broken, providing a solid relation between breaking 
of supersymmetry and that of U(1)h- 

2.4 Supersymmetry breaking and the low-energy spectrum 

To induce supersymmetry breaking VEV Fx, we need a linear term of X in the superpotential. 
This implies that U(1)h must be broken either explicitly or spontaneously. Here we simply 
parameterize the effect of U(1)h breaking in the X potential by a dimensionless chiral spurious 
parameter 77, which we assume to have the U(l) H charge of —2. The resulting physics does not 
depend much on the underlying origin of this breaking. 

The most general low-energy 4D interactions of X consistent with the (broken) U(1) H sym- 
metry is given by the following Kahler potential and superpotential: 

T^XtX- J-(XtX) 2 + ---, (32) 
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w ^ + ^ x+ ^ + wk> x3+ -' (33) 

where c is a constant term in the superpotential, needed to cancel the cosmological constant, 
and \j? x = Aur] A 2 . Here, again, we have used naive dimensional analysis to estimate the sizes of 
various coefficients (except for the c term), and omitted an 0(1) factor in each term. 5 Note that 
the terms in Eqs. (32, 33) arise from operators localized on the y = ttR brane, except for the c 
term which can have contributions from other sources as well. 

The scalar potential arising from Eqs. (32, 33) can be minimized in supergravity. Assuming 
that the coefficient of the (X^X) 2 /A 2 term in the Kahler potential is negative, the minimum 
of X is given by the competition between the X mass term arising from V ~ (fi x / A 2 )|X| 2 C 
\dW/dX\ 2 (d 2 K/dX^dX)~ 1 and the linear term V ~ — 2/j 2 x c(X + X^)/Mp X arising in supergrav- 
ity. The constant c is determined to cancel the vacuum energy V ~ \dW/dX\ 2 — 3|W| 2 /-^pi as 
c ~ /i^Mpi/V^. This, therefore, leads to the following supersymmetry breaking minimum 

2A 2 

W ^ ^Kfa « !0 14 GeV, F x ~ fx, (34) 

with the mass-squared for the X excitation given by m 2 x xs ix x j 'A 2 . Note that the X VEV, 
(X) ss 10 14 GeV, is smaller than the compactification scale, 1/nR 10 15 GeV, so that the 
4D analysis of the potential minimization is justified. In fact, with \xx much smaller than (X) 
to reproduce the weak scale superparticle masses (see Eqs. (28, 30, 31) and below), the only 
relevant terms in the potential minimization are the first two terms of Eqs. (32) and (33). 

The supersymmetry breaking of Eq. (34) can be transmitted to the MSSM gauginos and 
scalars by gauge mediation by coupling X to the messenger fields / and /: W — XX ff [22]. 
The minimum of X in Eq. (34) is not destabilized as long as the coupling A is sufficiently 
small, A 2 nj/167r 2 ^ (A/Mpi) 2 , where rif is the number of components for the messenger fields. 
We introduce the messenger fields in the bulk as hypermultiplets: {/, f c } and {f,f c }- The 
boundary conditions are given by 

/(5) = f D (Wt$ ® fL(l,2)%?\ (35) 
/ c (5*) = fU3\ltf ] ®m,2ti7 2 \ (36) 



and similarly for {/, f c }, leading to the zero modes from f D , f L , f D and f L . Here, we have chosen 
the messenger fields to be a pair of 5 + 5*, for simplicity, but they can in general be an arbitrary 
number of pairs of arbitrary SU (5) representations (as long as they do not make the standard 



5 The most general insertions of the spurious parameter rj allows us to write down the tree-level \x term in the 
superpotential, with fi w AtttjA w ^/A. This contribution is the same order as the one in Eq. (28); see Eq. (34). 
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model gauge couplings strong at or below ~ 1/-R). The messenger fields have interactions to X 
on the y = nR brane: 



S = Jd 4 xj* R dy5(y-nR) jd 2 6 ( Vd Xf D f D + i] L XfJ L ) + h.c, (37) 



10 

where the couplings rjo and t]l are of order An/M* (4it/M^ 2 ) from naive dimensional analysis 
if X is a y = nR brane (bulk) field. This determines the U(1)h charges of the f = Jd + II and 
/ = fo + jh fields such that the sum of the / and / charges is —2. (The f c and f c fields have 
the opposite charges to / and /, respectively.) 

The messenger multiplets in general have the bulk mass terms of the form of Eq. (12), Mf 
and Mf. The interactions of Eq. (37) then lead to the 4D superpotential 

W = X D Xf D f D + X L XfJ L , (38) 

where /l, fc and Jl represent the zero-mode chiral superfields, and 

\ D ps \ L m Airefej, (39) 

where e/,e/ are given by Eq. (17) with $ = /,/. The stability condition for the potential is 
\ 2 D L n//167r 2 ^ (A/Mpi) 2 pa 10~ 4 , which can be easily satisfied, for example, by taking Mf, Mj ^ 
1/irR, i.e., fn, Jl-, Jd and /l localized towards the y = brane. At the scale 

M mcss » A D , L (X) » (40) 

the messenger fields are integrated out, generating the gauge-mediated contributions to the 
MSSM gaugino and scalar masses [17, 18]: 



9a F * m 2 -2M Vr^'f 9 " Y 



(41) 



where a = 1, 2, 3 represents the standard model gauge group factors, g a are the standard model 
gauge couplings at M mcss , f — q,u, d, I, e, H u , H d , and C{ are the quadratic Casimir coefficients. 

The supersymmetry breaking parameters and the \x parameter in our theory receive contri- 
butions of Eqs. (28, 30, 31) generated at the scale M* and those of Eq. (41) generated at the 
scale M mess . The low-energy superparticle masses are then obtained by evolving the parameters 
of Eqs. (28, 30, 31) from M* to M mess , adding the contributions of Eq. (41) at M mess , and then 
evolving the resulting parameters from M mess down to the weak scale. Because of the wavefunc- 
tion suppression factors 6f j, which are exponentially sensitive to the bulk masses Mf j, the value 
of M mess can in general be anywhere between ps 100 TeV and O(0.1) (X) ^ 10 13 GeV. Here, 



12 



the upper bound comes from the stability condition on Xd,l, while the lower bound from the 
messenger stability. Note that since the gauge- mediated contributions of Eq. (41) have the size 

M _ ( m 2\i/2 „ F x ( 9 2 M Pl \ F x , , 

Ma ~ {m f } ~ X \J6^~A~ ) ~ X' (42) 

where (7 represents the standard model gauge couplings, they are comparable to the tree-level 
contributions to the Higgs-sector parameters of Eq. (28). 6 On the other hand, the flavor non- 
universal contributions of Eqs. (30, 31) are suppressed due to the e factors associated with 
the quark and lepton superfields (except possibly for the third generation). This therefore 
reproduces precisely the pattern for the low-energy supersymmetry breaking masses in flavorful 
super symmetry. 

The model also has other flavor violating contributions to the supersymmetry breaking pa- 
rameters, but they are all small. For example, loops of the higher dimensional gauge and 
messenger fields produce flavor violating scalar squared masses at 1/R, but they are of order 
N m _ css (g 2 /167T 2 ) 2 \Fx/A\ 2 pz ((X) / A) 2 m 2 j and thus small. They = brane Kahler potential opera- 
tors connecting the matter (and messenger) fields, e.g. S(y) J d 4 6 TjTjTlTi and 5{y) J d 4 6 T^Tjpf, 
also generate flavor violating scalar squared masses through loops of the matter (or messenger) 
fields. Using naive dimensional analysis to estimate the coefficients of the operators, we find that 
this contribution is at most of order \Fx / A\ 2 / (nRM*) 5 and negligible. Possible contributions 
from bulk higher dimension operators are also expected to be small based on similar dimensional 
arguments. Finally, y = brane localized kinetic terms, e.g. 5(y) Jd 4 9T^Tj, can introduce flavor 
violation by giving corrections of order 1/M*R pz 1/IQtt 2 to the kinetic terms of the low energy 
4D fields. After canonically normalizing the 4D fields, these corrections affect both the Yukawa 
couplings and the supersymmetry breaking parameters. Interestingly, however, this does not 
affect the mass insertion parameters used in section 3.1 at the order of magnitude level. In other 
words, we can always take the basis for the low energy 4D fields such that the Yukawa couplings 
and supersymmetry breaking masses are given by Eqs. (16, 30, 31) at M* even in the presence 
of the general brane kinetic terms. 7 Below, we assume that this basis is taken. 

Setting the size of the dominant contributions to the supersymmetry breaking and \x param- 
eters to be the weak scale, we obtain F x /A pz (100 GeV - 1 TeV) from Eq. (42). The value of 
F x is then determined as y/F^ pz (10 8 - 5 - 10 9 - 5 ) GeV using Eq. (27). This leads to the gravitino 
mass 

7713/2 ~ vSpi * (tU " 10) GeV ' (43) 



6 In contrast with the situation discussed in Rcf. [14], there is no reason in the present theory that the \i term 
must be suppressed compared with the gauge- mediated contributions. In fact, they are naturally expected to be 
comparable. 

7 In fact, this property persists even if the corrections to the 4D kinetic terms are of order unity. 
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implying that the gravitino is the lightest supersymmetric particle (LSP). Together with the 
flavor non-universal contributions of Eqs. (30, 31), this can lead to spectacular signatures at the 
LHC [2], some of which will be discussed in section 3.5. 

2.5 Neutrino masses, R parity, and dimension five proton decay 

The U(1)h charge assignment of Eqs. (24, 25) contains two free parameters qQ and qi- These 
parameters can be restricted by imposing various phenomenological requirements [15]. For ex- 
ample, if we require that dangerous dimension-five proton decay operators W ~ QiQjQkLi and 
UiUjD k Ei are prohibited by U (1)h, then we obtain the conditions SqQ+qL ^ and SqQ+qL ^ —4, 
respectively. Similarly, if we require that U(1)h forbids dimension-four _R-parity violating op- 
erators W ~ LiH u , QiDjL k , UiDjD k , LiLjE k and K ~ L\H d , we obtain q L ^ —1, q L ^ 1, 
QQ ^ -!» Ql ^ 1 and q L ^ 1. 

An interesting possibility arises if qi = 0. In this case we can have the following superpoten- 
tial on the y = nR brane: 

S = Jd 4 x £% 5{y - 7cR) J d 2 9 ^-XO N% N] + h.c, (44) 

which, together with the last term of Eq. (14), leads to 

W = ^>XN t N 3 + (yJijLiNjHv, (45) 

in the low-energy 4D theory. Using naive dimensional analysis, the couplings and (y u )ij are 
given by Kij ~ 4^6^^ and (y u )ij ~ 4n e^Ny The vacuum of Eq. (34) is not destabilized 
as long as ^ O(0. 1), which can be easily satisfied by taking to be somewhat small, i.e., 
by taking <: —1/nR. Small neutrino masses are then generated by the seesaw mechanism 
through the X VEV of Eq. (34). Note that the e Ni factors cancel out from the generated neutrino 
masses: 

Wij ~ 4neL * eL >i~xy' ^ 

It is interesting that with (X) 10 14 GeV, this is in the right ballpark to explain the experi- 
mental data on neutrino oscillations. 8 

It is not necessary to impose all the requirements above for the U(1)h charge assignment. 
For example, i?-parity violating operators can be forbidden simply by imposing matter parity 

8 The interactions of Eq. (45) also generate supersymmetry breaking masses of order (yl/lQir 2 )Fx / {X) for Li 
and H u through loops of Ni (A terms at one loop and non-holomorphic supersymmetry breaking masses at two 
loops [23]). This effect, however, is small for y v -C 1, compared with the contributions of Eqs. (28, 41). 
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in addition to U(1) H . Nevertheless, it is interesting that one can consider the U(1) H assignment 
that satisfies all these requirements. For example, one can adopt 

Qq = \ + 2n, q L = 0, (47) 

where n is an integer. The U(1)h symmetry is spontaneously broken by the VEV of X, but the 
charge assignment of Eq. (47) leaves a discrete Z 6 symmetry after the breaking. The product 
of Zq and U(l)y contains the (anomalous) Z% baryon number and (anomaly-free) Z 2 matter 
parity (R parity) as subgroups. This symmetry, therefore, strictly forbids the i?-parity violating 
operators, and the lightest supersymmetric particle is absolutely stable. 

In the rest of the paper, we assume that the LSP is absolutely stable (although it is not 
necessarily required by the model). This can be achieved either by choosing the U(l) H charges 
so that all the i?-parity violating operators are forbidden even after the U(1)h breaking, as is 
the case for Eq. (47), or simply by imposing matter (or R) parity. 

2.6 Origin of U(1)h breaking 

In section 2.4, we have simply parameterized the effect of (small) U{1)h breaking by a spurious 
parameter 77 <C 1. This breaking controls the size of the coefficient /jl x for the X linear term 
in the superpotential, and thus the size of supersymmetry breaking. There are a variety of 
possibilities for the origin of the required small breaking. For example, it may simply arise as a 
result of string theory dynamics at the cutoff scale M*. Here, we discuss two explicit examples 
for the origin of U(1)h breaking. The validity of the model as well as its basic phenomenological 
consequences discussed in section 3 have little dependence on this physics. 

The first possibility is that the U(1)h breaking effect arises from the mixed U(1)h anomaly 
with respect to the hidden sector gauge group. The scale fix then arises from dimensional 
transmutation associated with the hidden sector gauge group. This scenario can be implemented 
in our higher dimensional framework simply by promoting the model discussed in Refs. [15, 24] 
to higher dimensions. Specifically, we consider a supersymmetric S'?7(5)hid x SU (5) gauge theory 
on 5D flat spacetime, where the latter SU(5) factor is identified with the unified gauge group, 
whose gauge multiplet obeys the boundary conditions of Eqs. (1, 2). The Higgs and matter fields 
are singlet under S77(5)hid) and have the same SU(5) gauge quantum numbers and boundary 
conditions as in section 2.1. The location for the Higgs, matter and X fields, as well as their 
U(1) H charges, are also the same as before. 

The messenger fields {/, f c } and {/, f c } are also introduced in the bulk as before, with 
the interactions to the X field given by Eq. (37). Instead of introducing arbitrary explicit 
U(1)h breaking, however, here we assign the gauge quantum numbers (5*, 5) to / and f c , and 



15 



(5, 5*) to / and f c , where the numbers in parentheses represent the quantum numbers under 
SU (5) hid x SU(5). Below the compactification scale ss 1/irR, this reduces to the model discussed 
in [15, 24]. In particular, the required X linear term in the superpotential is generated: 

W cS = XA 2 hid X, (48) 

where Ahid is the dynamical scale of SU (5)hid, and we have taken A# m \l ~ A. Note that this 
superpotential is "exact," i.e., no higher order terms in X are generated. 

A virtue of the higher dimensional setup in the context of SU{b\iA x SU(5) is that the 
nontrivial wavefunction profiles of / and / needed to suppress \d,l (to satisfy the stability 
condition X 2 D L ^ 10~ 3 ) also suppress the superpotential coupling W = £ffH u H<i/A in the low- 
energy 4D theory, which can arise from the y = 7iR brane localized operator and leads to an 
unwanted large \i term unless ( ^ Xd,l- Using naive dimensional analysis, we find \d,l ~ C ~ 
47re/ej, so that we do not have a large ji term from the superpotential operator. 

Another possibility for the U(1)h breaking is that U(1)h is spontaneously broken. Since 
U(1)h has a mixed anomaly with respect to SU(3) C , this provides a solution to the strong CP 
problem [25]. We do not attempt here to construct a complete model of this kind. It is, however, 
straightforward to realize this possibility at the level of a non-linear sigma model, i.e. the axion 
field being realized nonlinearly. 

3 Phenomenology 

In this section we study phenomenology of the model presented in the previous section. We 
study constraints from flavor and CP violation and the variation of the superparticle spectrum 
allowed by these constraints. We find that there are a variety of possibilities for the next-to- 
lightest supersymmetric particle (NLSP), which decays into the LSP gravitino with the lifetime 
of O(10 2 - 10 6 sec). We also discuss proton decay, precision gauge coupling unification, and 
possible experimental signatures. 

3.1 Constraints from flavor violation and the variety of the spectrum 

Phenomenology of the model depends on the wavefunction profiles for the quark and lepton zero 
modes, which are controlled by the bulk masses for these fields. In the low-energy 4D theory, 
these affect the Yukawa matrices, Eq. (16), and the flavor violating contribution to the squark 
and slepton masses generated at M*, Eqs. (30, 31). This effect is parameterized by the factors 
e$ ($ = Q u Ui, Di, L i7 E h N>) in Eq. (17). 

The values for the e$ factors are restricted by requiring that the gross structure of the 
observed quark and lepton masses and mixings are reproduced by these factors. This, however, 
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still leaves some freedoms for the choice of the e$ factors. For example, scaling {eQ ; , ey v e^} — > 
{aeQ^ a~ 1 eu 1 , C(~ 1 ^d 1 } does not change the quark masses and mixings. Taking these freedoms 
into account, here we consider 

1 1 1 

eQ 2 ~F 5 a<?e, e f/ 2 ~ sT 3 **^, e» 2 ~ f ^^e, (49) 

~ e, e El w y'^a^ap e 2 , 

e L2 ^y~iaie, e E2 ~ y'^a^ap e, (50) 



with 



tan /3 ~ o;^ 6 1 , (51) 



where e = O(0. 1) and a q , ai and ap are numbers parameterizing the freedoms unfixed by the 
data of the quark and lepton masses and mixings. Note that the range of a. q ,i,p is restricted such 
that the e parameters, ^Q i ^j i) D i ,L i ,E i i do not exceed ~ 1; see Eq. (17). (The value of ap is also 
restricted so that tan/3 stays within the regime in which none of the Yukawa couplings blow 
up below the cutoff scale.) The pattern of Eqs. (49, 50) is a straightforward generalization of 
Eq. (18), and the resulting quark and lepton masses and mixings are still given by Eqs. (22, 23). 

The parameters a q , ai and ap, however, alter the size of the flavor violating contribution to 
the squark and slepton masses, and are thus constrained by low-energy flavor and CP violating 
processes. We use the mass insertion method [26] to derive constraints on these parameters. 
The experimental bounds on the mass insertion parameters can be found, e.g., in Ref. [27], and 
are summarized in Ref. [2]. In the quark sector, the most stringent bounds come from K-K, 
D-D and B-B mixings and sin 2(3, while in the lepton sector the most stringent one comes from 
the \i — > e7 process, giving 



^\Re(5 d 12 ) 2 LL/RR \ < (lO^-lO- 1 ), ^\Re(5 d 12 ) LL (5f 2 ) RR \ < 10" 3 , 

J\M&)l L/RR \ & (io~ 3 -io- 2 ), v /|im(^ 2 ) LL (^ 2 ) RR | <; io- 4 , 
m 2 )LL/ RR \ & (io-Mo- 1 ), |o&» ) L£ | = m 2 ) RR \ < (io- 3 -io- 2 ), 

\{^)ll,rr\ £ (0.1-1), \(5f 3 ) LL \ = \(5f 3 ) RR \ & 10- 2 , 



(52) 



(53) 



|(^ 2 ) iL |^(10- 4 -10- 3 ), (54) 
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where we have kept only the bounds relevant to our model. In deriving the above bounds, we have 
taken the gluino and squark masses to be the same order of magnitude with rriq ~ 500 GeV, 
and the same for the weak gaugino and slepton masses with ~ 200 GeV. For heavier 
superparticles, the bounds become weaker linearly with increasing superparticle masses, except 
for that on |(5i 2 )ll|, which scales quadratically with m ; ~. 

In order to compare our model with the above bounds, we need to obtain the structure of 
the squark and slepton mass matrices at low energies. We first consider the flavor universal 
contribution. It comes from two different sources. The first is gauge mediation, generated at 
the scale M mess , while the other is a U(l)y Fayet-Iliopoulos D-term piece, Tr(Y^ra^) 7^ 0, of the 
soft masses generated at M*, Eqs. (28, 30, 31). The sfermion masses at a low energy, can 
then be written as 



3 _„4 



a=l 



N 

1 ' mess 



1 + 1 - 




It 



SYffM f + f^_ 2 ^ + ^3 ln ^ > (55) 



5 16tt 2 V d y yoi\ yaf J A 2 m 

where (61,62,^3) = (33/5,1,-3) are the 321 beta-function coefficients, Yf represents hyper- 
charges in the normalization that Q has Yj = 1/6, and XH u ,H d ,Q 3 ,u 3 ,E 3 are the 0(1) factors in 
front of the corresponding soft masses generated at M*. (Here, we have kept only the leading 
terms in e.) As we will see in section 3.2, the U(1) Y -D-term piece can considerably affect the 
superparticle spectrum, leading to interesting phenomenology. 

The flavor violating elements of the sfermion mass matrices are renormalized among them- 
selves, and are also generated from the flavor universal piece through the Yukawa couplings. 
These effects, however, do not significantly modify the values of these elements in most of the 
parameter space. We therefore take the approximation that the flavor non-universal part of 
the sfermion masses is parameterized by Eqs. (30, 31) with Eqs. (49, 50) at low energies. 9 The 
chirality-preserving mass insertion parameters are then obtained by dividing these flavor violat- 
ing elements by the (average) diagonal elements in the super-CKM basis. 

With the low-energy mass parameters described above, one can study the constraints from 
flavor and CP violation. The scalar trilinear interactions in our model are generated only by 
renormalization group evolution, so that they are proportional to the corresponding Yukawa 
couplings with real proportionality constants, in the basis where the gaugino masses are real. 
They, therefore, do not contribute to flavor or CP violating processes. The constraints on the 
a parameters are then obtained from Eqs. (52 - 54). We find that for y = An and e = 0.05, all 



9 A possible contribution to m~ from loops of the right-handed neutrinos is also not important as long as 
{Vv)ij ^ O(l), which is the case for the e factor assignment of Eq. (50) with ai ss 0(1). 
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constraints from the quark sector are satisfied, while ji — > cy gives 



ol x <> 1.8, (56) 

with no further constraints on a q or ap. Taking e = 0.1, the constraints become stronger with 
both fj, — > e7 and K-K mixing, giving 

ai & 0.9, a p & 1.4. (57) 

These bounds are obtained for the superparticle mass scale of mj ~ 400 GeV, with Fx/M* ~ 
1 TeV. (This corresponds to m<j ~ 1.2 TeV, which is sufficient to avoid the LEP II bound on 
the physical Higgs boson mass.) While these bounds are rough ones, they show that there exists 
a consistent parameter region. For heavier superparticles, the bounds become weaker and the 
region expands. 



3.2 The NLSP 

As we have seen in section 2.4, the LSP is the gravitino with mass ~ (0.1 - 10) GeV. In order 
to study phenomenology, it is important to determine which particles can be the NLSP. Since 
the dominant contribution to the masses of most supersymmetric particles comes from gauge 
mediation, we first consider the spectrum without the corrections from tree-level pieces generated 
at M*. Since the masses are determined by the gauge charge, the lightest particles will be those 
neutral under SU(3)c and SU(2) L . Therefore, the lightest gaugino is a neutralino, Xi which is 
mostly bino, and the lightest sfermions are the right-handed sleptons. The mass of the bino at 
low energy is given by 

( \ at 9i(^r) Fx /_ q ^ 

m flN " ^mess , x y (58) 

while the mass of the sleptons can be derived from Eq. (55). From these two equations we see that 
with increasing N mess the sleptons become lighter than the bino, while increasing M mess makes the 
sleptons heavier because of renormalization group effects. Calculations show that for N mess = 1 
the bino is always the NLSP, while for larger N mcss the sleptons can be lighter. In the case of 
-^mess = 3 (5), for example, the sleptons are lighter than the bino for M mess ^ 10 10 (10 12 ) GeV. 

The bino mass in the present model is the same as in gauge mediation, but the slepton masses 
can deviate. As discussed in section 3.1, the sleptons receive the contribution from the U(l)y 
D-term, indicated by the second line of Eq. (55). This contribution is flavor universal so it does 
not affect the splitting among sleptons, but it affects the relation between the sleptons and the 
bino. The other correction to gauge mediation comes from the tree-level masses in Eqs. (30, 31). 
From Eqs. (49, 50), we see that these mass terms are e suppressed for most fields, but the effect 
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can be 0(1) for f R , and the unknown coefficient could even be negative as long as the sum of 
the tree-level and gauge mediated pieces bring the physical mass above direct detection bounds. 
This means that f R could lie anywhere in the spectrum of e R , jl R and B. 

The splitting between e R and j2 R is controlled almost entirely by the splitting at M* be- 
cause the renormalization group running is universal up to small effects from the muon Yukawa 
coupling. Phenomenology is governed by the splitting between mass eigenstates which is given 
by 

m L~ m L 0(0.01) /^X 2 /A/M P1 \ 2 m| 



2^4- n Lss uu o-oi J 

The splitting between light generation sfermions is much larger than in the usual gauge mediation 
scenario. It can be large enough that the heavier one can decay to the lighter by emission of an 
electron and a muon. 

There are corners of parameter space where the NLSP is not a right-handed slepton or bino. 
Since the contribution from the U(l)y -D-term in Eq. (55) has opposite signs for the left-handed 
and right-handed sleptons, it could invert the usual order between these two species. The lighter 
stop could also be the NLSP because, like f R , it has an 0(1) tree-level contribution to its mass. 
The stops also have a contribution from the large top Yukawa coupling, which decreases the 
masses through renormalization group evolution. While the tree-level piece is expected to be 
smaller than the SU(3)c gauge mediation piece, negative tree-level and Yukawa effects could 
combine to give a strongly interacting NLSP. We do not consider these exotic NLSPs in the rest 
of this paper because they require large cancellation between independent effects. 

3.3 Proton decay 

Dimension four proton decay in the present model can be forbidden by the U(1)h symmetry or 
matter parity. Dimension five proton decay caused by colored Higgsino exchange is also absent 
because of the form of the Higgsino mass matrix determined by higher dimensional spacetime 
symmetry [9] . Proton decay in the present model can thus arise only from dimension six operators 
and cutoff suppressed dimension five operators. 

As discussed in section 2.5, we can take the charge assignment of U(1)h such that the 
operators W ~ QiQjQkLi and UiUjDkEi are forbidden: SqQ+qL 7^ 0, —4. In this case, dimension 
five proton decay arises only from operators on the y = irR brane which involve the X VEV. 
The relevant interactions are W ~ X m QiQjQ k Li and X m U i UjD k Ei ) which can be written for 
3qQ + qL = —2m and 3gg + <?l = 2m — 4 (m G Z > 0), respectively. In the low-energy 4D 
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effective theory, these interactions lead to dimension five operators 

W a 4ne Q .eQ.e Qk e Ll -^-QiQ j Q k L l and 4 neu . eu . eDkeEi ——U i U j D k E l , (60) 

where the coefficients are evaluated using naive dimensional analysis, and we have used (X) m 
A 2 /Mpi. We find that the approximate sizes of these operators are obtained by replacing the 
colored Higgsino mass by 47rM^/A m_1 in the corresponding expressions in the minimal super- 
symmetric SU(5) grand unified theory. The resulting proton decay rate is thus much smaller 
than the current experimental bound for all the values of 3g<g + qi 7^ 0, —4. 

Dimension six operators are generated in the present model only through brane localized 
terms, since without them exchange of bulk gauge bosons does not transform a quark into a 
lepton or vice versa. (Note that different 321 multiplets arise from different SU(5) multiplets, 
see Eqs. (5 - 10).) The relevant terms are kinetic mixing operators K ~ T+T', T^T", F^F' and 
cutoff suppressed dimension six operators K ~ T^T'T^T", T^T'F'^F on the y = brane. Here, 
we have omitted factors involving the gauge multiplet needed to make operators gauge invariant, 
and the existence of Hermitian conjugates is implied. The kinetic mixing terms lead, through 
unified gauge boson exchange, to dimension six operators at low energies, whose coefficients 
have approximately the size obtained by replacing the unified gauge boson mass by 1 / nR in the 
corresponding minimal supersymmetric SU(5) expressions. For 1/irR as 10 15 GeV, this leads to 
a proton decay rate somewhat larger than the current experimental bound [28]. This implies that 
the compactification scale should be somewhat larger (by a factor of a few) or the coefficients 
of the original kinetic mixing operators should be suppressed (by an order of magnitude or so). 
This potential difficulty does not arise in weakly coupled models, an example of which will be 
discussed in section 5. The coefficients of low-energy dimension six operators arising from the 
cutoff suppressed operators are similar in size to those in the minimal supersymmetric SU(5) 
model, so that they do not lead to proton decay at a dangerous level. 

In summary, proton decay in the present model is caused by dimension six operators, orig- 
inating from terms on the y = brane. Since the wavefunction values for the first and second 
generation fields on this brane are typically of the same order, the proton can decay into fi- 
nal states containing /i + with a similar rate to those containing e + . This provides interesting 
signatures for future proton decay experiments. 

3.4 Precision gauge coupling unification 

Strongly coupled grand unification in higher dimensions allows a precise calculation for gauge 
coupling unification [12, 13]. Incalculable corrections arising from the cutoff scale physics are 
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suppressed, and the corrections from higher dimensional fields between the energy interval be- 
tween M* and 1/nR are precisely calculated. Here we study this issue in the model of section 2. 

We phrase the degree of the success of gauge coupling unification in terms of the prediction of 
a s (M z ) = g\{Mz)I^K obtained from gx^M-z)-, where 5^2,3 represent the standard model gauge 
couplings. In particular, we consider the deviation of the prediction in the present model, af D , 
from that obtained by assuming the exact unification in the MSSM, a^ GUT '°: 

5a s = af - af UT <° ~ -— a]A. (61) 

s s 27r s V 1 

Here, A parameterizes corrections from higher dimensional fields, which can be calculated within 
higher dimensional effective field theory. Using the result of Ref. [29], we find that in the present 
model 

3 15 
A = -- \n(nRM^) - 3 ln(e 0l eQ 2 eQ 3 ) + — ln^e^e^) 

9 9 6 

+- ln(e Dl e D2 e D3 ) - - ln(e Ll e i2 e L3 ) + - \n(e El e E2 e E . 3 ), (62) 

where we have used the approximation that the Higgs doublets are strictly localized to the 
y = tiR brane. (The term —(9/7) ln(e^e^ d ) should be added to the right-hand-side if the Higgs 
fields are delocalized.) Inserting Eqs. (49, 50) into this equation, we obtain 

3 135 45 45 

A — —— ln(7ri?M*) — In a q — — In ai + — In otp. (63) 

Considering that the logarithms are expected to be of order unity, we find that A is typically 
of O(10), with the sign depending on the values of a q ,i,/3- For typical superparticle spectra, 
including the one considered here, a good fit to the experimental values of <7i,2,3(Mz) is obtained 
for 

A cxp w 5 ± 0(1). (64) 
The expression in our model, Eq. (63), can easily accommodate this value. 

3.5 Collider signatures 

Phenomenology of the general flavorful supersymmetry scenario has been discussed in Ref. [2]. 
Here we summarize some of the basic features in the context of the present model. As we saw 
in section 3.2, this model has a large portion of parameter space where there is a charged NLSP 
which is stable for the purposes of collider studies. Unlike the conventional scenarios, the NLSP 
in flavorful supersymmetry could be a tr or a right-handed slepton of a different flavor. Heavy 
stable charged particles are relatively easy to see at colliders. By measuring their velocity and 
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momentum, their mass can be deduced. The mass of the charged NLSP be can measured to 
better than 1% at the LHC by measuing only a few hundred NLSPs with 0.6 < (3 < 0.91 [30]. 

Once the NLSP mass is known, it is possible to fully reconstruct events even in the hadronic 
environment of the LHC. Therefore we can determine the flavor content of the NLSP by taking 
its invariant mass with other leptons in the event. If the NLSP is found to be mostly selectron or 
smuon, this is definitive evidence for nontrivial flavor structure in the supersymmetry breaking 
sector, and possibly for flavorful supersymmetry. In addition, once we learn the dominant flavor 
of the NLSP, we can look for NLSP production in association with leptons of other flavors to 
measure the mixing angles of the NLSP. 

Because the lifetime of the NLSP is quite long, it can be studied in a cleaner environment. One 
proposal involves using the muon tracker to determine where in the surrounding rock an NLSP 
went, and extracting pieces of rock that likely contain NLSPs to study them elsewhere [31]. 
Another possibility is to build a large stopper detector outside of one of the main detectors 
which can stop the NLSPs and then measure the decay products [32]. This would allow precise 
measurements of the lifetime of the NLSP as well as the masses of the decay products. As pointed 
out in Ref. [2], a particularly distinct signature of flavorful supersymmetry is monochromatic 
electrons or muons in the decay of the NLSP, indicating a two body decay of a selectron or 
smuon. This is not a possibility in the conventional scenarios because the tr is the NLSP, and it 
decays to a r which further decays, so the many body decay causes the leptons to have a broad 
spectrum. Even if the NLSP is a f R , a stopper detector will allow us to look for rare decays 
into other flavors and precisely measure the flavor content of the NLSP. The stopper detector 
can also check to see if the LSP is the gravitino. From the kinematics, the mass of the LSP 
can be measured, which can then be tested against the supergravity prediction which relates the 
lifetime of the NLSP to the mass of the gravitino [33]. 

While the signatures are much more spectacular if there is a slepton NLSP, evidence for 
flavorful supersymmetry can still be found with a neutralino NLSP. One possibility is to look for 
direct slepton production from Drell-Yan processes and measure the spectrum through kinematic 
variables such as M^ 2 [34]. This is difficult because it requires high statistics and the Drell-Yan 
cross section falls rapidly with increasing slepton mass. Another possibility is to look for multiple 
edges in flavor-tagged dilepton invariant mass distributions as in Ref. [35]. This will allow us 
to find different flavors of sleptons if they are separated by more than a few GeV, which we 
would expect in flavorful supersymmetry. Finally, we could also study the spectrum of left- 
handed sleptons or even squarks to look for flavor non-universality. While these measurements 
are more difficult than those with stable sleptons, they could still provide information on the 
flavor structure of the supersymmetry breaking sector. 
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4 4D Realization — Model in Warped Space 



The model in section 2 has been formulated in flat space, but we can also consider a similar 
model in warped space, along the lines of Ref. [36]. An interesting feature of this model is that 
it allows for a 4D interpretation through the AdS/CFT correspondence, providing a picture of 
realizing flavorful supersymmetry in a 4D setup. 
Specifically, we take the metric 

ds 2 = e' 2ky r]^dx^dx u + dy 2 , (65) 

where k denotes the inverse curvature radius of the warped spacetime. The two branes are 
located at y = (the UV brane) and y = nR (the IR brane). The scales of these branes are 
chosen to be k pa 10 17 GeV and k' = k e~ wkR pa 10 16 GeV, respectively. The cutoff scale of the 
5D theory is taken to be M* pa 10 18 GeV. The gauge symmetry structure is as described in 
section 2; the bulk 577(5) symmetry is broken to 321 on the IR brane at y — ttR. The IR brane 
thus serves the role of breaking the unified symmetry. 

The configuration of the matter and Higgs fields is as described in section 2. The locations of 
these fields are controlled by the bulk masses, and the resulting Yukawa couplings are given by 
Eq. (16), where the e factors are given by Eq. (17) with M$ — > M$ — k/2. The analysis of U(1)h 
and supersymmetry breaking is as in sections 2.3 - 2.6. (Note that the cutoff scale on the IR 
brane is warped down to Ml = M^e~ 7TkR pa 10 17 GeV.) This leads to phenomenology discussed 
in sections 3.1, 3.2 and 3.5. Dimension four and five proton decay is negligible for the reasons 
described in section 3.3. Dimension six proton decay is also not dangerous as the unified gauge 
boson mass is now of order irk' pa 10 16 GeV. For gauge coupling unification, we can show, using 
the results of [37], that the threshold correction is still given by the formula Eq. (62). (Note that 
the contribution from the Higgs doublets to differential running shuts off above Ml, since these 
fields are localized on the IR brane.) The experimental values of the low-energy gauge couplings 
are thus successfully reproduced, as seen in section 3.4. 

The model described here has the following 4D interpretation through the AdS/CFT cor- 
respondence. At very high energies above k' pa 10 16 GeV, the theory is a 4D supersymmetric 
S77(5) x G gauge theory, where S77(5) is the unified gauge group and G some quasi-conformal 
gauge group. There are three generations of matter fields, 3 x (10 + 5*) of S77(5) (and possibly 
three right-handed neutrinos), but not the Higgs fields. There are also fields charged under G, 
some of which are charged under 577(5) as well. At the scale k' pa 10 16 GeV, the G sector 
deviates from the conformal fixed point, breaking the unified S77(5) symmetry to 321 by the 
gauge dynamics. It also produces the MSSM Higgs doublets and the supersymmetry breaking 
sector containing X as composite states. The effective theory below k' is thus the MSSM (and 
possibly three right-handed neutrinos) together with the supersymmetry breaking sector. 
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An important point is that the interaction strengths of the matter fields to the G sector are 
controlled by the dimensions of operators coupling matter to fields charged under G. In general, 
these dimensions are generation dependent. Moreover, since G is strongly interacting above k! , 
the anomalous dimensions for these operators can be large. As a result, the interaction strengths 
of matter to the G sector strongly vary between different generations, and since the Higgs 
doublets and X arise as composite states of G, the interactions of matter to these states show 
strong generation dependence. Since the origin of this generation dependence is common for the 
matter couplings to the Higgs fields (the Yukawa couplings) and to the X field (supersymmetry 
breaking couplings), the patterns of these two classes of couplings are correlated. The correlation 
is exactly the one given in Eqs. (16, 29), realizing flavorful supersymmetry. 

We have considered here a 4D theory in which the G sector is quasi-conformal and has a large 
't Hooft coupling above the dynamical scale, motivated by the warped space construction. The 
dynamics described above, however, are independent of these assumptions. The same dynamics 
can also be incorporated, in principle, in a purely 4D theory whose 't Hooft coupling is not 
necessarily large above k! . The quasi-conformal nature of the dynamics is also not essential. It 
will be interesting to construct an explicit example of purely 4D theory in which the G sector 
exhibits different renormalization group behavior, e.g. asymptotic freedom, above the dynamical 
scale A G « 10 16 GeV. 

5 Weakly Coupled (Non-Unified) Models 

In this section we present a non-unified model of flavorful supersymmetry in higher dimensions. 
Here we do not require that the theory is strongly coupled at the cutoff scale, nor that it possesses 
the U(1)h symmetry. Rather, we assume that certain operators are small at the cutoff scale due 
to ultraviolet physics. 

We consider a supersymmetric SU(3)c x SU(2)l x U(1)y gauge theory in 5D flat spacetime, 
compactified on an S 1 /Z2 orbifold: < y < irR. As in the model of section 2, the two 
Higgs doublets are localized towards the y = ttR brane, where supersymmetry is broken by 
the F-term VEV of a chiral superfield X. The matter fields are introduced in the bulk as 
hypermultiplets, whose zero modes Qi, Ui, Di, Li, Ei (and iVj) are identified with the MSSM 
matter fields. The wavefunction profiles of the zero modes are controlled by the bulk masses M$ 
($ = Qi, Ui, Di, Li, Ei, Ni), as seen in section 2.2. 

We do not require that the theory is strongly coupled at the cutoff scale M*, which is taken 
to be a factor of a few above 1/R. We then naturally expect that the operators located on 
branes have 0(1) coefficients in units of M*. This leads to the 4D Yukawa couplings of Eq. (15) 
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SU(3) C x SU(2) L x U(1) Y 




y = y = ttR 

Figure 2: A schematic depiction of the configuration for various fields. 

with 

(y u )ij ~ e Qi e Vj , (y d )ij « e Qi e Di , (y e )ij ~ ez,^, (y„)ij ~ e^e^, (66) 

at low energies, where e$ are given by Eq. (17). By choosing e$ and tan/? to be as given in 
Eqs. (49 - 51) with y — 1, this reproduces the gross structure of the observed quark and lepton 
masses and mixings, Eqs. (22, 23). 10 The configuration of the matter fields, as well as those of 
the Higgs and supersymmetry breaking fields, are depicted schematically in Fig. 2. 

The supersymmetry breaking parameters are generated through the interactions of the MSSM 
states to the X field on the y = irR brane. In the absence of the U(l) H symmetry, the super- 
potential operators W ~ XQiUjH u + XQiDjH d + XLiEjH d + XLiNjH u are not forbidden in 
general. These operators generate flavor non-universal left-right mixing terms for the squarks 
and sleptons that require relatively heavy superparticles to avoid the constraints from low-energy 
flavor and CP violating processes. Here we assume that these operators are somehow suppressed. 
We also assume that the direct \i term, W ~ H u H d , is absent. Note that these assumptions are 
technically natural because of the nonrenormalization theorem. The supersymmetry breaking 
parameters are then generated by the Kahler potential operators and C ~ / d 2 6 XW^Waa + h.c, 

10 Here we have assumed that the Majorana masses for N{ are on the y = ttR brane, but not on the y = 
brane. This can be realized, for example, by introducing the U(1)b-l symmetry broken on the y = irR brane. 
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where W" (a = 1, 2, 3) are the 321 gauge field strength superfields, giving 




Here, we have omitted 0(1) coefficients in each term, and (77$)^ ~ e^ i e^ j ($ = Q, U, D, L, E) are 
general complex 3x3 matrices. This gives a correlation between the Yukawa couplings Eq. (66), 
and the supersymmetry breaking parameters Eqs. (67 - 72), realizing flavorful supersymmetry. 
Note that because of the absence of a factor 4-7T in Eq. (66), the mass splittings between different 
generation sfermions in Eqs. (68, 69) can be larger than those in the strongly coupled case. 

The model has other flavor violating contributions to the supersymmetry breaking param- 
eters, but they can be controlled. For example, loops of the higher dimensional gauge fields 
produce flavor violating supersymmetry breaking masses at l/R, but they are not much larger 
than the tree-level masses in the parameter region considered, as long as the coefficients of the 
matter brane kinetic operators at y = are of order l/167r 2 M* or smaller. Note that this size 
of the coefficients is technically natural. The matter 4-point Kahler potential operators on the 
y = brane also give flavor violating contributions at loop level. They are, however, suppressed 
by a factor of l/(nRM*) 5 and negligible. Possible contributions from bulk higher dimension 
operators are also expected to be small. 

The compactification scale l/R in the present model is naturally of order the unification 
scale Mjj 10 16 GeV to preserve the successful supersymmetric prediction for the low-energy 
gauge couplings. In this case, the gaugino and sfermion masses are of order m as Fx/Mu while 
the gravitino mass is m 3 / 2 ~ F x /M Ph so that m 3 / 2 ~ (Mu/Mp^rh (1 - 10) GeV, leading 
to signatures discussed in section 3.5 with the NLSP being one of the right-handed sleptons. 
The compactification scale, however, can in principle take any value larger than of order a 
few TeV, in which case the gravitino may be (much) lighter. Note that the supersymmetry 
breaking parameters of Eqs. (67 - 72) are running parameters evaluated at the scale l/R. The 
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low-energy superparticle spectrum is obtained by evolving them down to the weak scale using 
renormalization group equations. 

Here we have presented a non-unified model of flavorful supersymmetry in 5D. It is, however, 
straightforward to make it a unified model, e.g., based on SU(5). We simply have to adopt 
the field content and boundary conditions of section 2.1 and follow the analysis above. To un- 
derstand gauge coupling unification, we need to assume that incalculable brane-localized gauge 
kinetic terms on the y = ttR brane are somehow suppressed (or universal), but dangerous proton 
decay can be easily suppressed, possibly by U(1)r symmetry [9]. 11 It is also straightforward to 
extend the model to higher dimensions. The only requirement is that the Higgs fields and the 
supersymmetry breaking field X are localized in the same place in the extra dimensions. 12 An 
advantage of such a setup is that we can suppress cutoff scale dimension-five proton decay oper- 
ators by localizing the Q,U,E and D, L fields in different subspaces in higher dimensions. These 
extensions allow us to realize flavorful supersymmetry in a wide variety of higher dimensional 
models, with varying spacetime dimensions, compact space geometries, and gauge groups. 

6 Conclusions 

In this paper we have presented explicit models of flavorful supersymmetry in higher dimensions. 
The basic idea is to localize the Higgs fields and the supersymmetry breaking field in the same 
location in the extra dimension(s). The interactions of matter fields to the Higgs fields (the 
Yukawa couplings) and to the supersymmetry breaking field (operators generating the super- 
symmetry breaking parameters) then receive the same suppression factors from the wavefunction 
profiles of the matter fields. This leads to a specific correlation between these two classes of in- 
teractions, realizing flavorful supersymmetry. The resulting phenomenology at future colliders is 
very rich, while stringent experimental constraints from the low-energy flavor and CP violating 
processes can all be satisfied. 

We have constructed a unified model of flavorful supersymmetry in 5D, in which the theory is 
strongly coupled at the cutoff scale. Supersymmetry breaking is mediated to the supersymmetric 
standard model sector by a combination of cutoff suppressed operators and gauge mediation. 
This model addresses various issues in supersymmetric unification. We have also presented a 
model in warped space, which allows us to obtain a picture of realizing flavorful supersymmetry 
in a 4D setup, through the AdS/CFT correspondence. Finally, we have discussed models which 

11 It is interesting to note that the 321 gaugino masses do not have to be unified at the unification scale even 
if the model is unified because the gaugino mass operators reside on the y = ttR brane, where the active gauge 
group is only 321 [9]. 

12 To be more precise, it is sufficient to require that the matter interactions to the Higgs and X fields are 
suppressed by common wavefunction factors, allowing the Higgs and X to propagate in different subspaces. 



28 



do not require that the theory is strongly coupled at the cutoff scale. This construction can 
be easily extended to a wide variety of higher dimensional theories, with varying spacetime 
geometries and gauge groups. 

It is interesting to note that the present setup is very generic in the context of a single extra 
dimension. If we want to explain the observed hierarchical structure of the Yukawa couplings 
by wavefunction overlaps between the matter and Higgs fields, the simplest way is to localize 
the Higgs fields to one of the branes and lighter generation matter more towards the other 
brane. Now, if the supersymmetry breaking field X is not localized to the same brane as the 
Higgs fields, interactions of lighter generation matter to X are not suppressed, leading to large 
flavor violating supersymmetry breaking masses. To avoid this problem, we need to localize X 
to the same brane as the Higgs fields (unless some other flavor universal mediation mechanism 
dominates). This gives the spectrum of flavorful supersymmetry. 

As the LHC will turn on this year, it is important to explore possible theoretical constructions 
and experimental signatures of supersymmetric theories. The models presented here provide an 
example in which the supersymmetry breaking spectrum can be a window into the physics of 
flavor in the standard model. If supersymmetry is discovered at the LHC, it will be interesting 
to see if the longstanding assumption of flavor universality holds, or if there is a richer flavor 
structure within the supersymmetry breaking sector. This structure could give us information 
about the physics of flavor which could lie at energy scales as high as the unification or Planck 
scale. 
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